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The severe geomagnetic effects of solar storms or coronal mass ejections (CMEs) are to a 
large degree determined by their propagation direction with respect to Earth. There is a lack 
of understanding of the processes that determine their non-radial propagation. Here we 
present a synthesis of data from seven different space missions of a fast CME, which 


originated in an active region near the disk centre and, hence, a significant geomagnetic 


impact was forecasted. However, the CME is demonstrated to be channelled during eruption 
into a direction + 37 10° (longitude) away from its source region, leading only to minimal 
geomagnetic effects. In situ observations near Earth and Mars confirm the channelled CME 
motion, and are consistent with an ellipse shape of the CME-driven shock provided by the 
new Ellipse Evolution model, presented here. The results enhance our understanding of CME 
propagation and shape, which can help to improve space weather forecasts. 
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weather and lead to massive disturbances in the solar wind 

plasma and magnetic fields through the inner heliosphere 
up to the heliospheric boundary!~?. During a CME, a mass of the 
order of 10!? kg is expelled from the Sun’s corona, its outermost 
atmospheric layer, with kinetic energies of around 10% J (ref. 1). 
CMEs play a pivotal role in solar and space physics. They are 
responsible for the strongest disturbances in the geophysical 
environment, potentially leading to power blackouts and satellite 
failures at Earth’. Increasingly, policy makers recognize CMEs as 
a serious natural hazard, and counter-measures for the protection 
of space and ground-based assets are implemented. 

A major requirement for producing reliable CME forecasts is 
to know their direction as accurately as possible as they propagate 
away from the Sun. Indeed, previous research>-? has found that 
various heliospheric structures may alter the CME trajectory to 
change its geomagnetic impact drastically. The strongest change 
in the propagation direction from its solar source position, which 
coincides with the flare for strong eruptions!®, has been mainly 
argued to occur within the first few solar radii where the magnetic 
forces acting on the CME are strongest”!!~}’, although possible 
changes of CME direction during interplanetary propagation 
have also been put forward”. Non-radial CME eruptions of up to 
25° in longitude were predicted® due to the channelling of the 
CME to the location of the streamer belt, where the coronal 
magnetic field strength has a minimum. Coronal holes, which are 
regions of high magnetic field strength from where the fast solar 
wind emanates, are able to deflect CMEs away from their source 
regions, depending on their area and distance to the CME 
source®!*+1>_ However, precise quantifications and the maximum 
possible amount of deflection remain unclear. Deflected 
motions from high to low solar latitudes have been described 
for prominences, which upon eruption often form part of a 
CME!617, 

A major obstacle for quantifying the CME deflection in 
heliospheric longitude is the lack of coronagraphs that can image 
CMEs from outside the ecliptic plane. Consequently, deflections 
in heliospheric longitude are much more difficult to analyse than 
those in latitude. To accurately quantify this process, it is thus 
necessary to study the complete chain of solar, coronagraphic, 
heliospheric and in situ observations. This is now possible with 
the multiple imaging and in situ spacecraft available, forming the 
Great Heliospheric Observatory. 

In this work, we discuss an event that emphasizes the pressing 
need for improved real-time predictions of the geomagnetic 
effects of CMEs. On 7 January 2014, a very fast CME erupted 
from a solar active region facing Earth. Fast CMEs that erupt 
from source regions close to the solar disk centre are usually 
expected to impact Earth!®!8, so many observers around the 
world predicted that this CME would be strongly geo-effective 
and yet, no geomagnetic storm followed. We demonstrate that 
this CME was strongly channelled into a non-radial direction by 
the effects of its locally surrounding magnetic field rather than by 
coronal holes (CHs). We also show how the CME evolved in the 
heliosphere up to its arrival at Earth and Mars, confirming the 
inferences from solar imaging. To this end, we introduce the 
Ellipse Evolution (EIEvo) model for studying the CME 
propagation, which lets us derive constraints on the global 
shape of the CME-driven shock. To explain what happened in 
this event, we take a tour from the solar observations into 
interplanetary space, providing a synthesis of observations from 
13 instruments on seven space missions. 


C oronal mass ejections (CMEs) are the ‘hurricanes’ of space 


Results 
Solar on-disk observations. On 7 January 2014 19 universal time 
(UT), a CME erupted from a source region at 12° south and 8° 


west (S12W08) of disk centre, accompanied by a plethora of 
phenomena such as a flare, coronal dimmings, a global coronal 
wave and post-eruption arcades. These are the classical on-disk 
signatures of an erupting CME)’. The flare, peaking at 18:32 UT, 
was of the highest class (X1.2). A very fast CME (projected speed 
of ~2,400kms~') was observed by the coronagraph 
instruments onboard the Solar and Heliospheric Observatory 
(SOHO) in real time. Alerts for a G3 class geomagnetic storm 
or higher (on a scale from G1 to G5) were sent out by various 
space weather prediction centres around the globe”! and were 
picked up by the media. 

Figure la shows the state of the solar corona on 7 January 2014 
19:30 UT, about 1h after the flare peak, imaged by the 
Solar Dynamics Observatory Atmospheric Imaging Assembly 
(SDO/AIA7“) at 193 A. Several active regions can be seen, with 
the largest one at the centre of the solar disk. Post-eruption 
arcades, magnetic loops filled with hot plasma and counterpart 
signatures of an erupting CME magnetic flux rope, are visible at 
the flare site, located at the southwest corner of the large active 
region. Two large CHs are visible in the northeastern quadrant. 
In Fig. 1b, the results from the automatic Coronal Pulse 
Identification and Tracking Algorithm”? visualize the location 
of the front of a global coronal wave, which is thought to be 
driven by the lateral expansion of the CME”*. The algorithm 
could successfully track the wave almost exclusively in the 
southwest quadrant of the solar disk. This is confirmed by a visual 
identification of the wave in running difference movies 
(Supplementary Movie 1) showing an asymmetric wave 
propagation. A similar pattern can be seen in the coronal 
dimming regions in Fig. 1c (derived? from SDO/AIA at 211 A), 
which emphasize the evacuation of the corona at the locations of 
the footpoints of the erupting flux rope!? (Supplementary 
Movie 2). The dimming appears earlier and is more cohesive in 
the southwest quadrant than in the northeast. Consequently, both 
the coronal wave and the dimming provide early suggestive 
evidence of the non-radial motion of the erupting CME to the 
southwest as seen from Earth. What could be the physical cause 
of this non-radial propagation? 


Influence of CHs and solar magnetic fields. In Fig. la, the 
locations of two large CHs in the northeastern hemisphere of the 
disk seem consistent with the hypothesis that the CHs acted as to 
deflect the CME into the opposite direction. We took a closer 
look at this hypothesis within the framework of the so-called 
Coronal Hole Influence Parameter (CHIP!*!>). The CHIP 
depends on the distance from the CH centre to the CME 
source region, the area of the CH as well as on the average 
magnetic field inside the CH. It can be considered as a parameter 
describing how strongly the CME is pushed away from the CH. 
The CHIP value for the CME event under study (see Methods) is 
F=0.9 + 0.2G, which is a factor of 3-5 below that necessary to 
deflect a CME originating close to the central meridian almost 
completely away from Earth!*!>, In particular, the distance of the 
CHs to the flare is comparably large!, which lets us dismiss the 
hypothesis that the CHs are mainly responsible for deflecting this 
CME. 

The explanation for the non-radial propagation may rather be 
found in the solar magnetic fields near the eruption. Figure 1d 
shows the line-of-sight component of the photospheric magnetic 
field (SDO/Helioseismic and Magnetic Imager (HMI)°). An 
active region (AR 11944) is right at disk centre, a few degrees east 
of the flare position. The positive polarity (white) sunspot had a 
particularly strong vertical magnetic field of ~3,000 G, which 
exceeds usual values?” by a factor of about 1.5. The flare is located 
in between this strong sunspot and the small negative polarity at 
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its southwest (Fig. 1d), and the erupting CME orientation (see 
next section) is consistent with the direction of the photospheric 
inversion line in between those two magnetic polarities. The 
negative polarity is almost surrounded by positive polarities so 
that a coronal magnetic null point and related separatrices are 
expected. The study of this magnetic topology is of primary 
importance to understand flare reconnection (for example, flare 
ribbon locations). A further study, outside the scope of the 
present paper, including magnetic field extrapolation and data- 
driven magnetohydrodynamic (MHD) simulations, would be 
needed to understand precisely the role of the AR magnetic field 
complexity on the early CME development. However, as this 
topology is local, it is unlikely to be important for the CME 
development on scales larger than the AR, which is the focus of 
the present paper. 

Next, looking at Fig. lef a potential field source surface?® 
model shows that the streamer belt of closed field lines is highly 
inclined with respect to the solar equator (typical of solar 
maximum), and runs from north to south right above the strong 
active region. The CME source region is not under the streamer, 
but close to an area of open flux further west (green field lines in 
Fig. le,f). This provides some evidence that the CME has erupted 
in the direction of least resistance in the solar global field’, 
consistent with results of numerical simulations!*!3. The solar 
observations thus imply that the strongly non-radial motion of 
this CME is due to a combination of two effects: (i) the strong 
nearby active region magnetic fields to the northeast, and (ii) the 
open coronal field to the west of the source. Both processes acted 
to channel the CME to the southwest of the solar disk, which was 
reflected in the asymmetries of the global coronal wave and 
dimmings. 


SDO/AIA 193 post-eruption 1930 UT 


Global EUV wave 1804-1819 UT 


Coronal evolution. We now take a look at multi-viewpoint 
coronagraph observations of the CME in Fig. 2a. We used two 
methods to estimate the CME propagation direction up to 30 
solar radii (Ro). The first is the Graduated Cylindrical 
Shell (GCS) model, by which a wire-grid of a tapered hollow 
tube is fitted onto coronagraph images”??? by manual variation 
of several parameters controlling its shape. The triple viewpoints 
from the STEREO-B (CORI1/2 (ref. 31)), SOHO (C2/C3) 
and STEREO-A (COR1/2) imagers constrain the results very 
well°°, At the time of the event, the two STEREO spacecraft were 
151° ahead (A) and 153° behind (B) in heliospheric longitude 
with respect to Earth, at distances of 0.96 AU (A) and 1.08 AU 
(B). The CME propagates to the east in STEREO-B, where the 
event is seen as backsided, which shows that the CME longitude 
must be greater than — 153° + 180° = 27° west of the Sun-Earth 
line. The GCS model was applied between 7 January 18:15 and 
19:30 UT, when the resulting model apex position was between 
2.1 and 18.5 Ro. The average three-dimensional speed of the 
CME apex was 2,565+250kms~!, derived from a linear fit to 
R(t), not far from the fastest speeds ever observed? 
(~3,000kms~!). A constant CME direction is consistent with 
the time evolution in the coronagraph images, which gives 
32+10° (west) and —25+5° (south, with quoted errors 
common for the method?™?). This means that already very 
close to the Sun, at 2.1 Ro, the final direction of the CME was 
attained. 

A second method was used to find the speed of the CME 
segment that propagates in the ecliptic plane. We applied a 
triangulation technique’? to the CME leading edge in SOHO/C2/ 
C3 and STEREO-A/COR2/HI] (ref. 34) observations. Our results 
are averages of two methods (Fixed-f and Harmonic Mean) 
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Figure 1 | Solar observations of the X1.2 flare and associated phenomena on 7 January 2014 18-20 UT. (a) Location of coronal holes and post eruption 
arcade in SDO/AIA 193 A. (b) Extreme ultraviolet (EUV) wave evolution between 18:04 and 18:19 UT derived with the Coronal Pulse Identification 

and Tracking Algorithm (CorPITA) algorithm. Colours indicate the position of the wave front at different times. (c) Final extent of the coronal dimming 
(SDO/AIA 211A). (d) SDO/HMI line-of-sight magnetic field, showing the position of the large active region and the flare. White (black) colours indicate 
positive (negative) magnetic field polarities. (e,f) Pre-eruption potential field source surface model of the solar global magnetic field, as seen from Earth (e) 
and 40° west of Earth (f). It depicts closed (white) and open field lines (pink negative polarity, green positive polarity). Solar east (west) is to the left (right) 


in all images. 
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Figure 2 | Graduate Cylindrical Shell model of the CME and 
interplanetary shock kinematics. (a) Fit of the Graduated Cylindrical Shell 
(GCS) model (blue grid) overlaid on multipoint coronagraph observations, 
from left to right: STEREO-B, SOHO, STEREO-A. Shown are results for 7 
January 2014 at 18:25 UT +1min, when the GCS apex was at 4.2 Ro. 
(b) Distance (top) and speed (bottom) of the CME shock in the ecliptic are 
shown as a function of time. Blue (orange) solid lines are the kinematics 
towards Earth (Mars) calculated with the ElEvo model, based on a DBM 
with y = 0.165 x 10-7 km7! and w=400kms~'. Blue and orange dashed 
lines indicate errors from a variation of y from 0.16 to 0.17 x 10-7km~1, 
which results from the uncertainty in tmars of +1h. Black triangles are the 
results of triangulation, with speeds and their errors deduced from a 
derivation of a spline fit on the distance measurements. The observed 
arrival times and speeds at Earth are indicated by black diamonds, and the 
arrival window at Mars with a black horizontal line. 


assuming a small and wide CME extent in the ecliptic along the 
line of sight? . From a linear fit to R(t) between 20 and 30 Ro, we 
find a speed of 2,124 +283 kms” !, slightly lower than the apex 
speed from the GCS method. The direction of the CME front 
above 20 Ro is W45 10°. This is further west from Earth by 
about 13° compared with the GCS results, which can be expected 
since parts of CMEs seen at different latitudes may travel in 
slightly different directions, because of the CME three- 
dimensional tube shape. This is reasonable, as Fig. 2a shows the 
CME to be oriented with a moderate inclination angle to the 
ecliptic. Further from the Sun, we tracked the CME ecliptic 
leading edge to about 25° elongation with the STEREO-A 
Heliospheric Imager (HI), and applied the Fixed-Phi-Fittin 
method**. This results in a speed of 2,131+210kms~ |, 
consistent with triangulation. We use this as the initial speed 
for further modelling of the shock evolution in the ecliptic plane. 
We also assume that the CME leading edge is representative of 
the position of the CME-driven shock, which has been confirmed 
with imaging of CMEs at large elongations from the Sun and their 
in situ observations”. For the CME direction, we take 45 + 10° 
west of Earth or 37+10° away from the source region in 
heliospheric longitude, exceeding expected values for non-radial 
CME eruption in the corona’. 


Interplanetary evolution. Figure 2b extends the R(t) and V(t) 
functions of the CME shock up to Mars (1.66 AU). The inter- 
planetary kinematics towards Earth and Mars are shown together 
with the arrival times at both planets, which will be discussed 
further below. We modelled the shock kinematics with the drag- 
based model (DBM*°), which analytically describes the 
deceleration of CMEs by using equations of aerodynamic drag. 
It gives similar performances for arrival time predictions of CMEs 
at Earth as numerical simulations*’. The DBM has two free, 
constant parameters: the drag parameter y and the background 
solar wind speed w. Parameter y contains information on the 
CME mass and size, the ambient solar wind density and 
the interaction between the CME and the solar wind*°. For the 
background solar wind, we use w= 400kms~ !, an average solar 
wind speed observed at Earth by the Wind spacecraft a few days 
around the time of the CME. This means that there are no high- 
speed solar wind streams west of Earth near the CME principal 
direction. A previous CME on 6 January 2014 was directed 
towards STEREO-A*®, and is not expected to influence the 
propagation of the 7 January CME towards Earth and Mars. Both 
inferences support the view that we can safely use a constant 
direction of motion and constant y and w. Parameter 
y may vary*® from 0.1 to 2x 10~7km~!. Because Mars is 
close to the apex of the ecliptic part of the CME shock (6° away), 
we can choose a value for y that makes R(t) consistent with the 
shock arrival time at Mars (tmars, see next section). This results in 
7 =0.165 + 0.005 x 10~7km~!, an expected value*® for a fast 
CME propagating into a slow, unstructured solar wind. The time 
tMars has an uncertainty of + 1h (see below), providing the error 
margin in y. We now constrain the interplanetary shock evolution 
further with a new model for the shape of the CME shock in 
combination with the available in situ observations. 

Multipoint in situ observations of interplanetary CME (ICME) 
signatures, at longitude differences*”” from 10° to the size of the 
CME shock of ~ 100° can give constraints on how a CME evolves 
in the interplanetary medium*!~*°. Few such events have been 
described in the literature, and thus the global shape and 
extension of CMEs are poorly known. In a fortunate coincidence 
on 7 January 2014, Mars was at a heliospheric longitude of 51° 
west of Earth. Judging from the CME principal direction, Mars 
and Earth should see the apex and flank of the CME, respectively. 
Figure 3a visualizes the position of the planets and the STEREO 
spacecraft in the ecliptic plane, together with the shock evolution 
up to its arrival at Mars (see Supplementary Movie 3). We use the 
new EIEvo model for describing the global shape of the 
interplanetary shock (see Methods section), which is based on 
statistics of single-point shock observations*”. Before we discuss 
the range of possible parameters and errors for this model, we 
take a look at the solar wind and planetary in situ data, serving as 
boundary conditions for ElEvo. 


Arrival at Earth. Figure 3b-g shows the solar wind magnetic 
field and bulk plasma parameters observed by the Wind space- 
craft at the L1 point near Earth. For a very fast CME from a 
source region near disk centre, such as the 7 January 2014 event, 
the expected in situ ICME signatures are a shock, followed by a 
sheath region of enhanced plasma density and temperature, 
and a magnetic flux rope with a size of the order of 0.1 AU in the 
radial direction*®. However, in situ observations are limited to 
those acquired along a spacecraft trajectory through the global 
CME structure, and for a glancing encounter, the flux rope is 
likely to miss the spacecraft??. On 9 January 2014 at 19:40 UT, 
an interplanetary shock hit Wind, causing a sudden jump 
in the solar wind speed from 390 to 465kms~ !. This time, which 
we label as tgarth, defines the CME shock arrival at Earth. 
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Interplanetary evolution of 2014 January 7 19 UT CME 
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Figure 3 | Ellipse Evolution (EIEvo) model for the CME shock in the heliosphere and near-Earth solar wind. (a) Heliospheric positions of various planets 
and spacecraft on 7 January 2014 at 19:00 UT. The shape of the CME shock given by EIEvo is plotted for different timesteps as indicated by colours. 
The model parameters (bottom left) are stated for the last timestep tmars = 10 January 22:30 UT. For the same time, the ‘shape constraint’ gives a window 
for the heliocentric distance of the ellipse segment along the Sun-Earth line, making the ellipse shape consistent with teartn: (b-g) Solar wind magnetic field 
and bulk proton parameters in near-Earth space (Wind SWE/MFI) for 9-11 January 2014. (b) Total magnetic field (black) and components Bx (magenta) 
and By (yellow); (e) magnetic field B, component (in Geocentric Solar Ecliptic coordinates); (d) proton bulk speed; (e) proton density; (f) proton 
temperature (black, expected temperature“. from the solar wind speed in blue) and (g) geomagnetic Dst index. The first vertical solid line from the left 
indicates the arrival of the shock, and the second vertical line delimits the end of the ICME sheath region, which does not seem to be followed by a 


magnetic ejecta. 


No other CME can explain the arrival of this shock. Figure 2b 
demonstrates that the arrival speed of the shock given by EIEvo at 
Earth indeed matches the in situ shock speed of 488kms~}, 
derived from the MHD Rankine-Hugoniot relations. The shock 
was weak, with a magnetosonic Mach number of 1.2. 
The orientation of the shock using the co-planarity theorem 
points in the radial direction away from the Sun and 26° 
northward, which is consistent with the main direction of the 
CME to the south of the ecliptic plane*. A sheath region 
followed, with an average speed of 417+ 20kms~ ! and elevated 
proton temperatures, extending up to 10 January 2014 06:00 UT. 
This region has a radial size of 0.104 AU, with a magnetic field of 
9.5+1.9nT, which is enhanced compared with average solar 
wind values*®. This magnetic field is relatively weak compared 
with the average field usually found in ICMEs, which contain flux 
ropes*®. 

In summary, the above indicates that this CME almost entirely 
missed Earth, because a shock-sheath pair is seen but not any type 
of magnetic ejecta. Figure 3g shows the corresponding Dis- 
turbance storm time (Dst) index, which is derived from a 
combination of equatorial ground station magnetometers around 
the world. The sheath region has a magnetic field and radial size, 
which is comparable to non-cloud ejecta*® at 1 AU, but due to the 
predominantly northward B, (Fig. 3c), Dst does not even reach 
levels typical of a minor geomagnetic storm (Dst< — 50 nT). The 
maximum of the Kp index was only 3, which is below the NOAA 
threshold for a G1 category geomagnetic storm. However, the 
CME lead to a major solar energetic particle event near Earth, 
which resulted in an S3 solar radiation storm on the NOAA scale 
from S1 to S5. 


It is also interesting to note that the particular orientation of 
the CME favours a miss of the flux rope at Earth too, because the 
east ‘leg’ is far below the ecliptic, as demonstrated by the SOHO 
image in Fig. 2a. Thus, in addition to the non-radial eruption, this 
orientation should have contributed to the false forecasts, because 
CMEs with a moderate to high inclination of their axes to the 
ecliptic plane have a small angular extent in the ecliptic*?, making 
it likely that the flux rope inside this CME has crossed the ecliptic 
to the west of Earth. As a consequence, the possibly strong 
southward magnetic fields of the CME flux rope have not 
impacted Earth’s magnetosphere at all. 


Arrival at Mars. Figure 4a shows an electron spectrogram by 
Mars Express (MEX) Electron Spectrometer instrument”. 
Enhanced fluxes of electrons, originating from both the planet’s 
ionosphere and the solar wind and indicated by red colours in 
Fig. 4a, are seen starting on late 10 January. They show a drastic 
change on 11 January, when the electrons became both more 
intense and reached higher energies. These data suggest that a 
CME is arriving on this day at Mars, because the CME sheath 
contains both denser and faster plasma than the surrounding 
slow solar wind and additional energy is added into the induced 
magnetosphere of Mars by the CME. The MEX data do not allow 
a more precise definition of the arrival time than 11 January + 1 
day, so we turn to data from the Martian surface by the Radiation 
experiment (RAD!) onboard Mars Science Laboratory’s 
Curiosity rover. RAD is able to observe Forbush decreases 
(FDs), a temporary decline in galactic cosmic ray intensity when 
an ICME passes a detector’, Figure 4b shows RAD count rates 
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Figure 4 | Observations indicating the CME arrival at Mars. 

(a) At Mars Express, the CME is observed by the Electron Spectrometer 
(ELS) as an increase in the electron magnetosheath and solar wind 
differential energy flux (colour coded) starting late on 10 January, with clear 
enhancements on early 11 January to late 12 January. The horizontal arrow 
bar delimits the interval of enhanced high-energy electrons (30-400 eV) 
on 11-12 January. (b) Counts of energetic particles per second by the RAD 
experiment on the surface of Mars onboard Mars Science Laboratory's 
Curiosity Rover. The high-energy solar energetic particle event stems from 
an eruption on 6 January. The CME of our study was launched from the Sun 
on 7 January and its shock hit Mars on 10 January 22:30 UT +1h, as 
indicated by the onset of a Forbush decrease of the cosmic ray flux. 


per second of energetic particles, which includes primary particles 
and secondary ionizing radiation created by solar and galactic 
cosmic rays in the Martian atmosphere and regolith. After a solar 
energetic particle event, related to a different CME*®, a return to 
normal counts is seen followed by a FD onset at tgp = 10 January 
2014 22:30 UT + 1h. What appears to be a single-step behaviour 
of the FD in Fig. 4b could indicate an arrival of only the ICME 
shock, and not the flux rope”, How is the time tpp related to the 
ICME shock arrival at Mars? At Earth, there is in general a 
relatively tight correspondence in timing’! between ICME 
arrivals at L1 and ground-based FD onsets, with the shock 
mostly arriving a few hours before the FD onset. Given that the 
corresponding FD onset at Earth (not shown) is within <2h of 
the shock arrival at the Wind spacecraft, the start of the FD at 
Mars can be reasonably assumed to coincide, within a 2h 
window, with the arrival at Mars of the presumably strong shock 
driven by this fast CME: tmars = tep- 


Global shape of the shock. As we now know the arrival times of 
the CME at Earth and Mars, we can go back to Fig. 3a, which 
shows the shape for the CME shock in the ecliptic plane given by 
EIEvo for equidistant timesteps. The model assumes that the 
CME shock propagates along a constant direction with a constant 
ellipse aspect ratio (a,), constant angular half width (4) and one 
main axis of the ellipse oriented along the radial direction from 
the Sun. It is also possible to calculate the speeds and arrival times 
at in situ locations along the ellipse front analytically (see 
Methods section). The evolution of the ellipse apex, which 
is the point of the ellipse farthest from the Sun along the 
CME central direction, is modelled with the DBM. From an 


optimization analysis (see Methods section) using the multipoint 
in situ arrival times, we find the ellipse aspect ratio to vary as 
a, = 1.4 + 0.4, for half widths of the ellipse ranging from 35 to 60° 
(under the condition that the ellipse always hits Earth) and for the 
central ellipse direction at W45 + 10°. Plotted in Fig. 3a is an 
ellipse with parameters a, = 1.43, = 50° and direction W45. The 
multipoint modelling of the CME shock thus implies an elliptical 
shape elongated perpendicular to the propagation direction’. 


Discussion 

We presented for the first time the full evidence for a very strong 
non-radial motion of a CME, with complete observations from 
the Sun and two planetary impacts. The CME from 7 January 
2014 19:00 UT almost entirely missed Earth despite its source 
being close to the centre of the solar disk. We attributed this to a 
non-radial CME propagation direction, which was attained very 
close to the Sun (<2.1 Ro), rather than to a deflection in 
interplanetary space. The observations do not show a 
‘deflection’, which implies a change in direction, but rather a 
‘channelled’ CME motion, which is non-radial starting already 
with its inception on the Sun. We found a surprisingly large 
magnitude of this channelling with respect to the source region 
on the Sun (+ 37+ 10° in heliospheric longitude), and a so far 
largely unrecognized process causing the channeling by nearby 
active region magnetic fields*? rather than CHs°. The 
observations emphasize the need to understand the interplay 
between the active region and global magnetic fields in order to 
better predict the direction of CMEs, and support previous 
studies, which derive altered CME trajectories from modelling the 
background coronal magnetic field”. 

We also showed suggestive evidence that the non-radial CME 
motion can be seen in extreme ultraviolet observations of the 
corona, which showed asymmetries in the global coronal waves 
and dimming regions with respect to the flare position. Because 
such asymmetries can also arise from the structure of the solar 
corona’, further research is needed on the possibility of 
diagnosing non-radial CME motions within extreme ultraviolet 
images. Finally, the arrivals of the CME-driven shock have been 
observed in situ by the Wind spacecraft near Earth and by the 
RAD experiment on Mars Science Laboratory’s Curiosity rover 
on the Martian surface. These observations, together with results 
provided by the new EIEvo model, show that these arrivals are 
consistent with the CME direction given by solar and 
coronagraph imaging. The in situ arrival times allowed to 
directly constrain the global shape of the CME-driven shock to an 
ellipse with aspect ratio of 1.4+0.4, with the ellipse elongated 
perpendicular to the direction of CME motion”. 

The enhanced understanding of non-radial CME propagation 
presented in our paper will be helpful for real-time CME 
forecasting!® in order to avoid false positive predictions, as it was 
the case for the event studied, that is, a CME that was predicted to 
impact Earth actually missed it. However, this means that a false 
negative forecast is also possible: a CME that is launched from a 
source region as much as 40° away in longitude from the Sun- 
Earth line may impact Earth centrally. It needs to be emphasized 
that non-radial CME motion in heliospheric longitude is very 
difficult to study because the images showing the CME radial 
distance from the Sun result from integrations along the line of 
sight. The upcoming Solar Orbiter mission®°, imaging for the first 
time the Sun and the heliosphere from outside of the ecliptic 
plane, will provide better insights into the dynamical processes 
responsible for non-radial CME eruptions. 

In summary, the presented observations demonstrate the high 
value of many different instruments in spatially distant locations 
to study solar storms. We were able to draw a consistent picture 
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Figure 5 | Derivation of the geometry of a self-similar expanding ellipse in the EIEvo model. (a) A CME leading edge, a shock or the front of a 
flux rope, described as an ellipse, propagates away from the Sun in the ecliptic or solar equatorial plane with constant angular width and aspect ratio. 
(b) Geometry for deriving the speed of any point along the front of the CME leading edge. 
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Figure 6 | Optimization of the EIEvo model shape with multipoint in situ 
observations of a CME shock arrival. The average residual distances (D) 
between the ellipse and the heliocentric distances of Earth and Mars, at the 
observed in situ shock arrival times, are plotted as function of the ellipse 
aspect ratio a,, for half widths of 46-60° heliospheric longitude as indicated 
by the colours. The CME direction is set to W45, the mean value derived 
from triangulation. 


of the evolution of a CME from its inception on the Sun, during 
which the event under study experienced a strongly non-radial 
motion, to the impacts at planets and spacecraft. These 
fundamental results should help to improve the reliability of 
real-time forecasts of space weather. 


Methods 


Calculating the CH influence parameter. The CHIP is given as 


p= (34) X 


with <B> the average magnetic field inside the CH, A the CH area corrected for 
projection, d the distance from the CH centre to the eruption site and ep the unit 
vector along this direction. This vector defines the direction in which an erupting 
CME will be deflected due to the presence of the coronal hole, under the condition 
that F is sufficiently large. The values for A and the barycenter position (the centre 
of the CH region weighted by pixel intensity) were provided by the Heliophysics 
Event Knowledge base”, obtained with the SPoCA algorithm*® used on SDO/AIA 
193 A images”. We calculated the distance from the CH barycentre to the eruption 
site (S12W08) with the great circle distance on the solar sphere. The areas of the 
northern (southern) CHs are 1.2+0.04 x 10! km? (0.3 + 0.01 x 10!! km2), the 
distances to the source are 7.0 x 10°km (4.5 x 10°km) and the average magnetic 
fields are 1.38 (3.76) G (measured from SDO/HMI’®). As there are two CHs 
present, the CHIP is treated as a vector, which is summed up for all CHs present on 
the disk!4. 
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The resulting CHIP from both CHs at the location of the source region at 
S12W08 is F= 0.9 + 0.2 G, with the uncertainties arising from the calculation of the 
area with different thresholding techniques and an uncertainty in the 
determination of the magnetic field. The direction ep arising from the CHIP 
analysis is towards a position angle (PA) of about 230° (PA is measured from solar 
north at 0° to east—left side in a solar image—at 90°, south at 180° and west at 
270°), which also seems at first glance consistent with the CME propagation 
direction to the southwest, but the low CHIP!+!5 means that the combined coronal 
holes are not sufficiently close, have a large enough area or strong magnetic field to 
explain the non-radial propagation of the CME. 


Derivation of the EIEvo model for the evolution of CMEs. To model the shape of 
the shock driven by the CME in the interplanetary medium, we introduce a new 
method that describes the shock as an ellipse in the ecliptic plane, which we call the 
ElEvo model. It allows the extension of one-dimensional models for heliospheric 
CME propagation*®°, which provide the distance-time R(t) and speed-time V(t) 
functions of the CME front, into two-dimensional models of the evolution of CME 
boundaries in the ecliptic plane with only a few lines of code. ElEvo can thus be 
used to visualize the shape of a CME shock between the Sun to a given planet or 
in situ observing spacecraft using fully analytical formulas. Further, the speed and 
arrival time of any point along the ellipse can be calculated analytically. 

EIEvo is an extension to a model describing CME boundaries as self-similar 
expanding (SSE) circles®*!, Whereas SSE has been designed to derive CME 
parameters from observations of CMEs at large angles from the Sun with a 
heliospheric imager**®!, it can also be used to propagate a CME into the 
heliosphere and calculate expected planetary arrivals and speeds®°. The initial 
speed, direction and width of a CME, which are known from coronagraph 
observations, can be used as initial conditions. The advantage of ElEvo is that the 
shape of an ellipse is more flexible than an SSE circle, and thus, better suited for 
consistent modelling with multipoint in situ observations because the aspect ratio 
of the ellipse is a free parameter. 

The assumptions of the ElEvo model are: (i) the angular width in heliospheric 
longitude of the CME boundary remains the same for all times, (ii) one principal 
axis of the ellipse is oriented along to the propagation direction, and (iii) the ellipse 
aspect ratio and (iv) direction are both constants. However, it is possible to change 
the aspect ratio and direction in a code as a function of time, but this is not 
implemented in the current study. For describing the interplanetary deceleration of 
the CME, we use the DBM°®, which describes the kinematics for the distance R(t) 
and the speed V(t) of the CME as function of time (t). It has two free parameters: 
the drag parameter y (on the order of 0.1 to 2 x 10~7km~!), which describes the 
amount of drag exerted by the solar wind on—strictly speaking—the CME flux 
rope, and w, which is an average of the background solar wind speed. By choosing 
low values of y it is possible to describe the shock propagation up to 1 AU with 
DBM, so it can be used to calculate CME shock arrival times and speeds. In 
summary, the EIEvo model in combination with the DBM has four free parameters: 
(i) the inverse aspect ratio f, (ii) the ellipse half angular width 4, (iii) the drag 
parameter y and (iv) the background solar wind speed w. We now derive the 
equations necessary to code the geometry of the EIEvo model to visualize the SSE 
ellipse as it propagates away from the Sun as well as the speed for each point along 
the ellipse front. 


Visualizing a self-similar propagating ellipse. Figure 5a shows the geometry 

of an ellipse under the assumptions described above. The R(t) of the ellipse apex 
(the point of the ellipse farthest from the Sun along the ellipse central direction) is 
given by DBM”. In this section, we derive equations for the ellipse semi-major axis 
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a and semi-minor axis b as a function of R(t), the inverse aspect ratio f and the half 
width 2. We use f= b/a rather than a, = a/b because it simplifies the following 
calculation. The equations 


f =b/a 
BHA (2) 


0 = arctan(b’/a’ tanp), 


follow from the definition of f and the definition of angle f, which is the angle 
between the semi-major axis a and the normal to the tangent at point T (Fig. 5a). 
The location of T is the point of tangency on the ellipse for a line originating at the 
Sun. It can be easily seen that f = 2 by checking the sum of the angles of the small 
orange triangle in Fig. 5a in relation to a larger triangle (not highlighted) 
containing the angles 4 and 1. The polar angle of the ellipse 0 is given by a 
relationship from general ellipse geometry between f and 0. It is important to 
emphasize that we further construct the ellipse based on this particular value of 0, 
for which a line with distance r connects the ellipse center C to point T. 
Combining equation (2) gives a relationship for the polar angle 0 based on 
known parameters, 


0 = arctan(f? tan 4). (3) 
From the law of sines on the large orange-shaded triangle in Fig. 5a we derive: 
sind sing 
san 3 4 
r R(t) — b (4) 


Angle « follows from the angles of the orange-shaded triangle, and distance r from 
the definition of an ellipse in polar coordinates: 


a= 90deg +0 — 2, 


3 5 (5) 
r = ab/\/(bcos0)” + (æ sinb)“. 
The last equation can be rewritten with the definition of f as 
o= 2 — 1)cos?0 + 1, 
ae (6 


r=b/o. 
Introducing ~ from equation (5) and the last equation for r into equation (4) then 
eliminates the unknowns (a, r) and expresses b in function of known variables: 


_ R(t) sind 

~ cos(A— 0) + asin2’ 
a=b/f, 0) 
c= R(t) —b. 


Equations (7) are the final description of the ellipse parameters. The minor axis b 
of the ellipse depends on all known variables (R(t), f, à) through 0 and œ, from 
equations (3) and (6). The major axis a then simply follows from the definition of f 
in equations (2). The heliocentric distance of the centre of the ellipse is parameter c 
(Fig. 5b), closing the model equations necessary for visualizing the ellipse. 


Calculation of speeds along the ellipse front. For comparison to in situ 
observations, which give parameters such as the speed and arrival time of the CME 
shock with very good accuracy“$, one needs to know the speed of any point along 
the ellipse front as a function of the ellipse parameters. This problem has been 
solved analytically for the circular SSE geometryć?, and we introduce here the 
corresponding analytic solution for ellipses. 

Figure 5b demonstrates the geometry, with A being the known angle between 
the CME central direction and for example Earth, which could also be any other 
planet or spacecraft in the solar wind. The direction of the apex with 
respect to a coordinate system including the Sun and Earth depends on 
different methods for CMEs observed with coronagraphs?™? or heliospheric 
imagers®), 

We introduce a coordinate system centred on the ellipse (Fig. 5b), with 
coordinate X being perpendicular to the CME propagation direction and Y 
orthogonal to X. Here, c is the vector from the Sun to the ellipse centre, d connects 
the Sun to the front edge of the ellipse in the direction of Earth (that is, d stops at 
the ellipse boundary and does not connect the Sun to the planet), and r connects 
the centre to the end point of d on the ellipse: 


c= (0,c), 
8 
d=c+r. (8) 
The problem consists in finding the norm of d as a function of A. There are two 
crossings of d with the ellipse, one at the rear and one at the front (Fig. 5b), which 
will form the two solutions of the problem. The coordinates of r can be expressed 
with the projections of the vector d-c in the X/Y coordinates: 


Then, rx and r, can be introduced into the definition of an ellipse in cartesian 


coordinates, 
Tx\ 2 ry\? 
= 2) a1. 
(2) +@) 
This results in the following expression, which was simplified with the definition of 
f from equations (2): 


(10) 


af sin? A + (dcosA —c)?= b°. (11) 
This quadratic equation gives two analytic solutions of the front and rear crossings 
of d with the ellipse as a function of the parameters b, c from equation (7), f from 
equation (2) and A: 


_ ccosA + \/(b? — c?)f? sin? A + b? cos? A 13 
f> sin? A + cos? A ` (12) 

The solution with the positive sign in front of the root is the ‘front’ solution (d,) 
and the one with the negative sign the ‘rear’ solution (d2). The speed V4 (t) of the 
ellipse at the position defined by the angle A is derived from the self-similar 
expansion of the ellipse, which implies a constant half width 1. The assumption of 
self-similar expansion means that the shape of the ejection must not change in 
time, so the ratio between speeds and distances for all points along the ellipse must 
be constant™: 


dı k 


V(t). (13) 


Further, the time when d;(t) is equal to the heliocentric distance of the planet or 
the in situ observer gives the arrival time of the ellipse at the in situ location, and 
the speed Va (t) at this time the arrival speed. 


Analysis of the aspect ratio and width of CME shocks. In this section, 

we discuss how to find optimal solutions for the EIEvo shape when multipoint 
observations of the ICME arrival are available. We first create a shock apex 
kinematic R(t) with the DBM with parameters y = 0.165 x 107 7km~! and 
w=400kms~!, which yields a DBM arrival time at Mars consistent with the 
observed arrival time tmars: From the R(t) apex, we calculated with equation (12) 
for a range of half widths from 45° <4 < 60° the parameter dı parth for the longitude 
of Earth, at Earth arrival time tgarth. This range for values of / is chosen because the 
half shock extension is thought to be around 50° in heliospheric longitude*®. From 
the in situ observations, we know that the shock has impacted Earth. Thus, for an 
ellipse apex direction of W45 + 10°, 2 must be larger than 45° + 10°, or the shock 
would not reach Earth. This also means that the half width (2) of the shock is > 35° 
in the ecliptic, consistent with the value of ~ 50°. 

Because the ellipse shape needs to be consistent with both Earth and Mars 
arrival times, we repeated the same procedure for Mars and calculated the average 
of the residual distances D between the ellipse front (d, for the corresponding 
values of A for Earth and Mars) and the heliocentric distances of both planets 
(dearth and dmars) at the respective observed arrival times tparth and tars: 


Adgarth = [dearth ga di garth |, 


Admars = |dars = di Mars|; (14) 


D= (Adgarth F Admars)/2. 


Figure 6 shows parameter D as function of a, varied from 0.5 to 2.5 and for 
different 2 indicated by colours. This plot was made with the average CME 
direction of W45. For each A, an optimal solution for a, exists where D has a 
minimum. At the minimum, the ellipse impacts Earth and Mars at the observed 
arrival times. For two in situ spacecraft that observe an ICME arrival, there are 
optimal solutions for pairs of A and a, that match the two given in situ arrival times 
for the CME shock if either å or a, is known or assumed. We note in passing that 
for CME events where three or more in situ arrival times would be available, there 
would be higher residuals, and a stronger constraint might be derived. The 
important new result is that 2 and a, are not independent of each other, and the 
result of the optimization procedure is a range of 1.69 >a,> 1.23 for half widths of 
45°<1<60° when keeping the direction constant at W45. For a half width of 50°, 
the optimal a, = 1.43 (used in Fig. 3a). In general, as the width increases the aspect 
ratio must become smaller to be consistent with the in situ arrival times. 
To fully include the errors from the direction determination, we also experimented 
with varying the shock apex direction by + 10° around W45, a typical error for 
CME directions by the triangulation method?’. For the W55 case, 55°<1<60°, 
the optimal aspect ratio is in the range 1.84>a,> 1.58. For the W35 case, 
35°<1<60°, the aspect ratio is 1.51>a,>0.99. Thus, including the errors from 
both the direction and varying the half width within reasonable values leads to a 
considerable possible variation in <a,> = 1.4 + 0.4. However, <a,> indicates the 
global shock shape to be a slightly elongated ellipse, being only close to a circular 
shape for a very extreme choice of parameters. 

To better visualize this optimization process, we illustrated in Fig. 3a a distance 


r=d k window named ‘shape constraint’ along the Sun-Earth line, where the ellipse at 
dx = dsin(A), Mars arrival time has to pass through in order to be consistent with the shock 
d, = dcos(A), (9) arrival time at Earth. After impacting the Wind spacecraft, the shock travelled 26h 
A 50 min + 1h until its apex hit Mars. We assume that during this time the shock 
: travelled with a speed of 488 km s~ !, which is close to the slow solar wind speed so 
ty = dye. we do not expect much deceleration. With such a speed, the shock was at tyars at a 
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distance of 0.315 + 0.015 AU further away from Earth, along the Sun-Earth line. 
This distance has an error (indicated on the figure by small horizontal lines on 
Fig. 3a) due to the uncertainty of + 1h in fmars: One can see in Fig. 3a that 

the outermost (red) ellipse indeed crosses the ‘shape constraint’ window, 

which means that the implementation of ElEvo is consistent with the observed 
multipoint arrival times. 
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